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ABSTRACT 

We study the statistical properties of the cosmological 21 cm signal from both the in- 
tergalactic medium (IGM) and minihalos, using a reionization simulation that includes 
a self-consistent treatment of minihalo photoevaporation. We consider two models for 
minihalo formation and three typical thermal states of the IGM - heating purely by 
ionization, heating from both ionizing and Lya photons, and a maximal "strong heat- 
ing" model. Wc find that the signal from the IGM is almost always dominant over 
that from minihalos. In our calculation, the differential brightness temperature, 5Ti,, 
of minihalos is never larger than 2 mK. Although there are indeed some differences in 
the signals from the minihalos and from the IGM, even with the planned generation of 
radio telescopes it will be unfeasible to detect them. However, minihalos significantly 
affect the ionization state of the IGM and the corresponding 21 cm flux. 



1 INTRODUCTION 



The 21 cm signal from neutral hydrogen has the potential to 
provide a direct and independent probe of the cosmological 
gas distribution and its thermal and ionization states before 
and during reionization. Indirectly, such observations can be 
used to study the formation and evolution of early luminous 
sources, their feedback on the intergalactic medium (IGM), 
and their impact on the formation of subsequent generations 
of structures. As the 21 cm line is associated with the hyper- 
fine transition in the ground state of neutral hydrogen, it de- 
pends on the distribution of atoms in these two energy levels, 
which is quantified by the spin temperature Ts . Since all the 
atoms are immersed in a sea of CMB radiation, in principle 
they can absorb or emit photons with the frequency corre- 
sponding to the energy splitting between the two hyperfine 
levels, vo — 1.4204X 10"^ MHz, possibly changing the shape of 
the CMB spectrum itself. If only CMB photons are present, 
absorption and emission equilibrate in a short time, and no 
line will be visible. However, there are other two mechanisms 
that can change the spin states of neutral hydrogen atoms 
- collisions and Lyg excitati on (Wouthuysen-Field effect, 
IWouthuvsenlll952l : lFieldlll958l ). which render the 21 cm line 
visible in emission {Ts > Tcmb) or absorption (Ts < Tcmb) 
again st the CMB background (e.g. iMadau. Meiksin fc ReesI 
1 19971 ). 

The 21 cm signal from the high redshi ft universe 
has b een studied in relation to minihalos (llliev et al.l 
l2002l : iFurlanetto fc Obi 120061: IShapiro et all |2006| ) and 



the IGM (e.g . 



iTozzi et al.l l2000li ICiardi fc Madai 
Furlc 

Mellema et all l2006l : 



He et al.l |2004 : IFurlanetto. Zaldarriaga fc Hernau"is3 

u 



Santos et al 



200 i1 ). 



20031: 



2004 



formation of radiation sources, the IGM is cold and neutral 
and the 21 cm line is expected to be observed in absorption 
against the CMB. Observations of the 21 cm signal at these 
very high redshifts would provide a wealth of information 
about the initial density fluctuations and the presence of 
a running spectral index or deviations fr om gaussianity 
jLoeb fc Zaldarriagall2004l: IAU. Bharadwai fc Pandavi,2005i: 



Pillepich. Porciani fc Matarrese _ 20071: Lewis fc Challinorl 



20071 : ICoorav. Li fc Melchiorrill2008l : iMao et al.ll2008h . Once 
the ionization sources turn on, they affect the 21 cm 
signal both by reducing the amount of neutral hydrogen 
and by increasing the temperature of the IGM though 
photoheating, Lya photon scatteri ng and, if the sources are 
hard e noug h, X-ra y heating (Madau, Meiksin fc Ree^ 
19971 : lOhl l200l| : |Chen fc Miralda-Escudc 2004 



Chen fc Miralda-Escudc 2008; Chuzhov fc Shapiro! 



Furlanetto fc Pritchard ,2006.: .Chuzhov fc Shapird 



Prior to the 



Ciardi fc Salvaterral l2007l: ISemelin. Combes fc Baelj I2OO7 
Ripamonti. Mapelli fc Zaroubil |2008| ). Eventually, such 
heating renders the IGM visible in emission. 

A 21 cm signal is also expected from minihalos (MHs), 
collapsed structures with masses larger than the Jeans mass 
and smaller than the mass corresponding to a virial temper- 
ature of ~ 10*K. In such halos, the gas cannot cool efficiently 
unless a large fraction of molecular hydrogen is present. 
This, together with their fragility wit h respect to feedback 
effects (see e.g. lCiardi fc Ferrarall200"5h makes their contribu- 
tion to star formation and photon production small at best. 
Nevertheless, MHs are expected to form copiously through- 
out cosmic history and their detection would be an obser- 
vational breakthrough. It has been suggested th at a possi- 
ble way of detecting MHs is via the 21 cm fine (H iev et al.l 
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\200i : iFurlanetto fc Ohll2006l : IShapiro et allbood ) . Because 
of their high density and temperature, MHs result in a signal 
that is expected to be different from that of the IGM. On 
one hand, collisions are frequent in such halos, coupling Ts to 
the gas temperature; on the other hand, the very dense gas 
can trap ionization fronts, protecting their central regions 
from being ionized. This makes MHs a possibly significant 
or even dominant source of 21 cm emission at some stage of 
cosmic history. 

As MHs are expected to be found preferentially in high 
density regions, they are more likely to absorb UV radia- 
tion produced by nearby sources and get progressively pho- 
toevaporated (Shapir o et al. 1997 : Shapiro fc Raga 2000l: 



Haiinan et al.i ,2001.: .Barkana fc Loebl 12002 



Shapiro et al 



Whalen et al 



lUev et al.ll2005l : lAhn fc Shapirol 120071 : 

As a consequence, they can deplete the number of 
ionizing photons and delay IGM reionization. In addition 
to photoionization, MH formation a nd evolution is affected 
by a variety of feedback effects (see ICiardi fc Ferraral |2005| 
for a complete review), who se relative importan ce has not 
been clearly established yet. I Ciardi et al ] (|2006l . C06 here- 
after) performed the first simulations of reionization with a 
self-consistent treatment of MH photoevaporation, finding 
that, depending on the strength of the feedback effects, the 
presence of MHs could delay the end of reionization by as 
much as Az « 4. In this paper, we use the results of these 
simulations to study the 21 cm signal from both the IGM 
and minihalos. 

The layout of the paper is as follows. In Section[2]we give 
a short description of the C06 simulations and the method 
used to calculate the 21 cm signal from the IGM and the 
MHs. In Section [3] we present our results and discuss their 
uncertainties. In Section |4] we give our conclusions. 



2 METHOD 

In this Section we describe the method used to derive the 
21 cm signal from the IGM and the MHs. 



2.1 Simulations of reionization with minihalo 
photoevaporation 

Here we give a short description of the simulation run in 
C06. Readers can find more details in the original paper 
and references therein. 

The underlying N-body simulation was run in a cu- 
bic box of comoving length L = 20/i~^Mpc, which 
was extracted and resimulated from a much larger re- 
gion (jYoshida. Sheth fc Diaferiol [SooJ). Galaxy formation 
was followed with a semi-analytic model (|Kauffmann et ahl 
1 19991 ) and radiative transfer was calculated with the code 
CRASH (ICiardi et al.l [ioOl: Maselli, Fcrr ara fc Ciardil l2003l : 
iMaselli. Ciardi fc Kanekaili2008 ). with the assumption that 
the gas follows the dark matter. This distribution was 
then mapped on a 12 8^ grid, using the Triangular Shaped 
Cloud interpolation l|Hocknev fc EastwoodI llQSlI ) to get 
the density in e a ch ce ll. Throughout the simulation, the 
lEfstathiou et al.1 l|l992l ) transfer function was used and 
the cosmology was taken to be (fim, f^A, ^6, /i, o"8, J^) = 
(0.3,0.7,0.04,0.7,0.9,1), where D.^, D,a, and ilt are the 
mean matter, vacuum, and baryonic densities in units of 



the critical density, h is the Hubble constant in units of 100 
kms~^ Mpc~^, o"8 is the rms amplitude of matter fluctu- 
ations on the 8 Mpc h~^ scale, and n is the slope of the 
primordial power spectrum. 

Unlike the original simulations 

by ICiardi. Stoehr fc White! (|2003D and 

ICiardi. Ferrara fc Whitel l|2003l ). in C06 sub-grid physics 
was included to take into account the effects of MHs on 
reionization. Generally speaking, the box size of a reioniza- 
tion simulation should be large enough to be representative 
of the whole universe. This requirement makes it imprac- 
tical to resolve individual minihalos (whose mass may be 
as small as 10^ Mq). For this reason, C06 used a "sub-grid 
correction" to include their effects on reionization, which 
are mainly due to two issues: 

1) Gas accretion onto MHs, which reduces the density 
of the IGM; 

2) MHs absorption of photons that would otherwise ion- 
ize the IGM. Because of their high gas density and recom- 
bination rate, MHs can consume many more photons than 
the IGM, relative to their overall mass fraction. 

To take into account the above effects, at each step of 
the simulation in C06, the following "sub-grid correction" 
was applied: 

• In each cell, the total mass of MHs was c alculated using 
the Extended Press-Schechter (EPS) theory (iLacev fc Colg 
1 19931 : IMo fc Whit j 1 19961 : ISheth fc Tormenl 120021 1 H. and the 
IGM density was corrected accordingly. 

• The number of ionizing photons absorbed by MHs in 
each cell was taken to be: 



(1) 



where tmh is the optical depth of MHs, which is proportional 
to the mass fraction of MHs /coIi.mh, and A'^, is the number 
of photons injected into this cell. 

• At each step, the fraction of minihalos photoevaporated 
was computed as: 



MH.Gva — 



A^7,MH 



(2) 



where nu is the hydrogen number density, Al is the size of 
one cell, while ^ is the average number of photons consumed 
by minihalos per total atom (the expressions can be found 
in the original paper). 

As mentioned in the introduction, the detailed effect of feed- 
back on MH formation and evolution is still unknown. For 
this reason COG have considered two extreme models in the 
simulations. The first is the extreme suppression (ES) model, 
in which once a cell has been crossed by an ionizing photon, 
no new halos are allowed to form within it unless it recom- 
bines to become completely neutral. The second is the refor- 
mation model (RE), in which MHs continue to form undis- 
turbed by feedback effects. The real situation is expected to 
lie between these two extreme cases. 



^ A slightly different prescription is used in cells hosting sources, 
but their impact is minor. 
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2.2 21 cm signal from the IGM 

The differential briRht ness temperature 5Tb o f the IGM 
is calculated foUowinK [Ciardi fc Madad (120031') (see also 



Furlanetto fc Loeb|[2002l : lMellema et al.ll2006l : IShapiro et all 
20061 '): 



STk = 0.016T 1 



^c^ 



where 



T= -(1+5)(1-/co11,mh)(1-2:) 



0.02 



(3) 



1 + z 
10 



0.3 



1/2 



(4) 

with X ionization fraction. This formula is accurate only in 
the optically thin limit, which is appropriate for the IGM. 
Ts is the spin temperature, which determines the distribu- 
tion of neutral hydrogen atoms in the two hyperfine levels. 
The value of Ts is a weighted mean between the CMB tem- 
perature, TcMB, and the gas kinetic temperature, Tk, and 
depends on the effect of collisions and Lya pumping. In this 
work we use 



TcMB + {ya,eff + yc)Tk 
1 + yc,eff + Vc 



(5) 



where yc 
ing H-H 



IS 

H-e 



the coUisional coupling effic iency ( i nclud- 
and H-p collision s; see iNussej I2OO5I : 



iKuhlen. Madau fc Montgomervll2006h . and ya,eff is an ef- 
fective Lyg coupling effici ency term that can be written as 
l|Chuzhov fc Shapirdliooel ): 



= exp [-0.3(1 + zy/^T,-^/''] ( 1 + ^ 



(6) 



where Va is the sta,ndard Lya coupling efficiency 
IjMadau. Meiksin fc Reeslll997l ). The IGM temperature in 
the absence of heating mechanisms has been calculated us- 
ing RECFAST (|Seager. Sasselov fc Scottll 19991 . 2000). 

Typically, the IGM can be found in three thermal states, 
which define the value of the differential brightness temper- 
ature, as listed below: 

1) No heating. In this state, which occurs before and 
during the early stages of structure evolution, the IGM cools 
adiabatically. As long as the density of neutral hydrogen is 
high enough, collisions between different atoms couple Ts 
to Tk, which is usually lower than Tcmb at that stage, so 
that the 21 cm signal is in absorption. However, following the 
expansion of the universe, the density decreases as oc (1-f 2)'', 
so that coUisional coupling becomes ineffective below z ~ 20. 
In this case, the spin temperature is only slightly smaller 
than TcMB . 

2) Lya and UV photon heating. Once the first sources 
of radiation turn on, their UV photons start to ion- 
ize and heat the surrounding gas. "While their impact is 
limited to the regions in the immediate vicinity of the 
sou rces, the Lya photons (both continuum and injected; 
see IChen fc Miralda-Escudej 12004 IChen fc Miralda-EscudI 
I2OO8I ) can travel much further from their production sites 
because of the longer mean free path compared to UV pho- 
tons, and quickly build up a background that can influence 
the IGM as a whole. The Lya photons have a double effect 
on the spin temperature. First, even a small Lya intensity 
(fti 10~'^°erg cm~^s~^Hz~^sr~^) can make ya in eq. ([6} much 
larger than 1 and bring Ts very close to Tk , so that the spin 
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Figure 1. Evolution of CMB temperature (long-dashed line), 
IGM temperature in the absence of heating mechanisms (solid), 
with LyQ heating (short-dashed), spin temperature without 
(dashed-dotted line) and with (dotted) Lya coupling. Lyo pho- 
tons are included only from the start of the simulation, i.e. from 
z = 18.5. 



temperature traces the state of the IGM. This effect is typ- 
ically referred to as Lya coupling or the Wouthuysen-Field 
effect. 

Secondly, Lya photons can heat or cool the 
IGM slightly by resonance sc attering, as has been 
studied by several authors ( Madau. Meiksin fc Rees 



1997 



Chen fc Miralda-Escude 



2004 



Rvbicki 



20061: iFurlanetto fc PritchardI I2OO6I ). Recen t works 
jChen fc Miralda-Escudd |2004 iRvbickil I2OO6I : 
IFurlanetto fc PritchardI bood : iMeiksin 2006) find lower 



values than originally estimated ( Madau, Meiksin fc Reej 
I1997D . Nevertheless, it may still ha ve some impact on 
the thermal evolution of the IGM (|Ciardi fc Salvaterral 
l2007f ). In this paper we calculate this eff e ct wit h the Lya 
background derived in ICiardi fc Madaul ^Of), which is 
self-consistent with the simulation of COG. Figure [1] shows 
the temperature evolution of a homogeneous cosmic density 
field. The dashed-dotted line is the spin temperature 
without any heating mechanisms, while the dotted line is 
the spin temperature used in our Lya heating case. Both 
estimates exclude the photoionization heating. "With the 
adopted Lya flux background, the temperature of the IGM 
is never higher than Tcmb before the end of reionization. 
Thus, as Lya heating is not extremely efficient, we expect 
a 21 cm signal in absorption. 

3) Strong heating. Finally, when the Lya background 
becomes strong enough or additional sources of heat- 
ing are present, we expect the temperature T^ to be- 
come much larger than Tcmb. Such high temperatures 
could be reached in the presen ce of X-rays, that can 
heat the IGM above 1000 K (Ch en fc Miral da-Escudj 
2004l:IChen fc Miralda-Escude 2008; F urlanetto fc Oh 200l; 



Furlanetto. Oh fc Briggsl I2OOQ ). Although there are stiU 
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some uncertainties on the extent of the X-ray background, 
it seems plausible to have Ts 2> Tcmb at later times dur- 
ing cosmic evolution. In this situation, the 21 cm signal 
is always in emission |Zaldarriaga. Furlanetto fc Hernguistl 



l2004l : ISantos. Coorav fc Knoxll200d ). 



In addition to X-rays, Lya and UV photons, shocks 
due to structure formation can also heat the IGM 
dFurlanetto fc Locb 2004 ; Minia ti et aLl l2004': Shapi ro etal] 
I2OO6I ) . Shock hea ting may be significant if the gas is still cold 
at high redsh i ft jFurlanetto. Oh fc Brigd [200^ ) . However, 
I Shapiro et alj (|2006l ) showed that even taking into account 
shock heating, the 21 cm signal from MHs is still larger than 
that from the IGM at 2; < 20. In fact, shock heating usually 
is efficient around large halos or in high density regions. The 
same regions however are also likely to be ionized by nearby 
sources. This may suppress the effect of heating by shocks. 
Also, C06 did not consider shock heating because they as- 
sumed the gas to trace the dark matter distribution. So this 
contribution is not considered further in the present work. 

The details about the transition between these three 
configurations are still not clear, mainly because of our per- 
sistent ignorance on the details associated with the heating 
mechanism. Here, we estimate the 21 cm signal in three heat- 
ing case: pure UV heating, Lya and UV heating and strong 
heating. The UV heating is treated as did in the simulation 
by C06. 



2.3 21 cm signal from MHs 

In principle, the differential brightness temperature of MHs 
can be derived by replacing T in eq. (j4]| with: 



Tmh = — (1+5)/co11,mh(1— a;) 



0.02 



1 + z 
10 



0.3 



1/2 



(7) 



which is proportional to the mass fraction of MHs. 

However, this formula is valid only in the optically-thin 
case, while the optical depth of an individual minihalo is 
usually large. In fact, for a halo with mass 10^ M©, the typ- 
ical value of the optical depth al ong a line-of-sigh t across 
the center of the halo is « 0.5 l|lliev et al.ll2002l ). and it 
drops as the mass increases because of the high Ts and the 
low a b sorption co e fficien t of high mass halos (see llliev et al] 
|2002| ). i lev et al. I (|2002l ) give the more appropriate result 
of 21 cm flux and STi, for an optically-thick case, assum- 
ing that the halo is a Truncated Iso thermal Sphere (TIS; 
IShapiro. Iliev fc Ragalll999l : llliev fc S hapiro 2001). In this 
paper we make use of their results. 

Because of their high density and temperature (com- 
pared with Tcmb ) , individual MHs ar e expected to alway s 
produce a 21 cm signal in e mission (IShapiro et al.l |2006| ). 
with a large 5Tb (in Fig. 1 of llliev et al.ll2002l . the typical 
value is several hundred mK at frequency vo = 1.4204 x 
10^ MHz, and even larger after integration along the line 
profile). However, it is extremely difficult to resolve indi- 
vidual halos, because their angular size is usually less than 
1 . For this reason it is customary to estimate, rather than 
the signal from individual halos, the beam-averag ed effective 
differ ential brightness temperature as (eq. (6) in llliev et al.l 
|2002| ): 

- 1^ ^'^^ff^T.^oA^dM, (8) 



where dn/dM is the mass function of minihalos and H{z) is 
the Hubble parameter. Here we take the Jeans mass as the 
lower mass limit, which is 



5.7 X lO-" 





2 




5 


0.15 




0.02 


V 10 ) 



(9) 



as in lShapiro et al.l (j2006h . The upper mass limit is Afmax — 
2.8 X 10^(1 + z)-^/^A/q, which corresponds to T^ir = W K 
in our cosmology. A = ivr^js is the geometric cross section of 
a halo, where rris is t he truncated radius in the TIS model 
(jlliev fc ShapirdbOOll ): Aty^aiz) = [(27r/x)^''2i/ocr„/c]/(l -f 2) 
is the redshifted effective line width and ^T^,^,, is the dif- 
ferential brightness temperature of a single halo at the 
frequency vq (t he detailed expressions can be found in 
llliev et al.ll2002l ). 

While dn/dM can be derived from a Press-Schechter 
approach, in this case, we need also to account for photo- 
evaporation of the MHs. While the simulations described 
in Section [5] provide the total MHs mass fraction in each 
cell at any given time, we have no information on the dis- 
tribution o f mass amongst t he MHs. We thus proceed as 
follows. In llliev et al.l l|2005l ) it is shown that the mass 
evolution of an individual halo exposed to a flux Fo — 
F/{10^'^s"V[47r(lMpc)^]}, with an initial mass Mo at red- 
shift zo, is a function M = M{Mo, zo, z, Fq). The function 
can be inverted to derive Mo = Mo{M, zo, z, Fo), i.e. from 
the mass of a halo we can get the mass of its progenitor. 
Then, at any redshift and in any given cell, the minihalo 
mass function, including photoevaporation, can be written 



dn 

1m 



iM,z) 



d^n 
, dM'^dz 



j(Mo,z ) X 



dn 
dM, 



{Mo,z, 



dM. 
1m 

.dMo 
dM 



^{m'o,M,z' ,z,Fo)dz' 

(Mo,M, Zcr,Z,Fo). 



(10) 



Here, Zcr is the redshift at which a photon packet crosses 
the cefl for the first time and m'q = Mo{M, z , z, Fo). The 
first term of the equation takes into account the evolution 
of the mass function due to all the halos forming between 
Zcr and z, while the second term describes the evolution 
in the mass function of halos already in place at Zcr- Both 
d^n/ dModz{Mo, z) and dn/dMo{Mo, z) are calculated with 
the EPS theory, following C06. The equation above is valid 
in the reformation case, while in the extreme suppression 
model, once a cell has been crossed by ionizing photons, no 
new halo can form. In this case, the first term of the equation 
should be set to zero. 

Throu ghout the paper, we use the mass photoevapora- 
tion law of llliev et al] (|2005| ). which is: 



(11) 



M{t) = Mo 1 



In this equation, C and B are parameters determined by 
the spectr um of the sources. In our case C = 1.07 and 
B = 2.5 (|lliev et allboOSh . The evaporation time tev — 
teviFo, Mo, Zcr) IS a fuuctiou of Fo, Mo, Zcr and of the source 
spectrum. Figure [5] shows the mass evolution of individual 
halos with different initial masses. Mo, at zo = 17.6, ex- 
posed to various constant fluxes of Fo= 0.01, 1, and 1000. 
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Figure 2. Redshift evolution of the remaining mass of halos with 
different initial mass, under the effect of a photoevaporating flux 
with Fq = 0.01, 1, 1000 (from top to bottom in each line group). 
The initial (at zg = 17.6) mass of the halos, Mg, is as labelled. 



Figure 3. The evolution of the mass function under photoevap- 
oration for the extreme suppression model (thin lines) and the 
reformation model (thick lines). From top to bottom, the lines 
correspond to redshift 17.6 (when the evaporation begins), 17.4, 
17.2 and 17.0 respectively. Here we assume a constant flux Fa = 1. 



While low-mass halos disappear on very short time scales 
(«14 Myrs if Fo = 1 and Mq = 10* Mq), larger mass halos 
retain most of their mass. For example, with Fo = 1, a. halo 
with Mo — 10* M0 will have a mass M = 6 X lO'^ Mq after 
126 Myrs, i.e. has undergone a mass loss of only 40%. 

In principle, if the flux evolution in each cell is known, it 
is possible to re-construct the mass evaporation process for 
halos in each mass bin according to eq. (|ll|l . In C06 , the flux 
is determin ed by sources in the simulation box (Ci ardi et al.l 
l200ll .[ 2006l l . However, as it is computationally unfeasible to 
record the value of Fq in each cell at each timestep of the sim- 
ulation, we only recorded the total evaporated mass fraction 
of minihalos in each cell. We then adopt an effective ioniza- 
tion flux which gives identical evaporated mass fraction as 
Fo, to calculate the mass evolution of halos with different 
mass. Two methods are used to get the effective flux and 
implement eq. HI}. In the first, in each cell and at each 
time step of the simulation, , we use an effective ioniza- 
tion flux, F^^^ff, which is assumed to be constant between 
Zcr and . As we know the total mass of MHs in the cell, 
Mtot ~ X], use the following equation: 



= ^Mo,, 1^1 



(12) 



where the sum is performed over all the MHs in the cell. 
The above equation is solved to derive the only unknown 
quantity FQ ,,ff, which is then used in eq. (llOp to calculate 
the re-distributed MH mass function and STf At the follow- 
ing time step, the same process is repeated to derive a new 
effective flux. 

In the second method, we also choose an effective flux 
for each cell. Differently from the first method, here the ef- 
fective flux evolves step by step, rather than being constant 



from Zcr to z*. At each time step, z'', we calculate: 



(13) 



where: 



am; = Mo,: 



Ctev{FQ^^j, J,, Moj, Zcr) 



(14) 



is the mass of minihalo j evaporated between two steps. 
From the above equations we derive the only unknown quan- 
tity, Fgcff and then we use it in eq. (fT(7|l to calculate the 
re-distributed MH mass function and STt- Also, the same 
process is repeated to derive a new effective flux at the fol- 
lowing time step. We have used both methods and found 
that the results are consistent with each other, as in the 
vast majority of cells at any time the difference in the STb 
obtained with the two methods is less than 5%. Thus, in the 
following we will always use the flrst one. 

Figure |3] instead shows the mass function at different 
redshifts with the assumption of a constant ionization flux 
Fq = 1, as calculated in eq. pO|l . The thin lines refer to 
the extreme suppression model, while the thick lines are the 
reformation model. As expected, the mass function in the 
reformation model is always larger than that in the extreme 
suppression model. In the former, at high mass end, the 
formation rate of halos is higher than the evaporation rate, 
so the mass function increases with time. 

Using the above mass functions, we can now calculate 
the expected differential brightness temperature. Figure |3] 
shows the mean (volume-averaged) STi, evolution for the two 
MH formation models, which has been calculated both with 
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Figure 4. Evolution of the volu me-averaged STf, o f minihalos 
from the simulations described in lCiardi et al.l ||2006|) . The solid 
(dashed) lines refer to an optically thick (thin) case. The thin 
(thick) lines refer to the extreme suppression (reformation) model 
for minihalo evolution. 



the optically-thin formula given in eq. ((3| with the assump- 
tion that Ts S> TcMB (dashed lines) and for the optically- 
thick case i n eq. (El) (solid l ines), with STt^^g taken from 
Figure 1 of llliev etall l|200^ . Although the STt for an in- 
dividual halo is very high compared with the optically-thin 
case, the beam-averaged STt is very small (with a typical 
value of « 1 — 2 mK). This is because the halos usually have a 
small geometric cross section A, and thus the associated flux 
is not very significant. Figure!?] also shows that the optically 
thin assumption with Ts S> Tcmb is a good approximation 
to calculate the mean 21 cm signal from minihalos. STt in 
the extreme suppression model drops dramatically following 
the evaporation of the minihalos. In the reformation model 
instead, it increases continuously until late times. Thus, the 
presence or absence of MH reformation can result in totally 
different pictures of the evolution of the 21 cm signal from 
MHs. 



3 RESULTS AND ANALYSIS 

In this Section we present our results in terms of mean dif- 
ferential brightness temperature (as described in Sees. 2.2 
and 2.3) and fluctuations of brightness temperature, and we 
discuss their uncertainties. As a visual example, in Figure [5] 
we show maps of the STt in the pure UV heating case and 
reformation model. Form top to bottom, the three rows cor- 
respond to the IGM, MHs, and the total. The map shows 
that at the early and late stages of reionization, MHs con- 
tribute more to the total signal compared with the IGM, 
while at intermediate redshifts the contribution from par- 
tially ionized IGM increases. 



3.1 Mean differential brightness temperature 

The contribution from MHs to the mean differential bright- 
ness temperature is calculated as described in the previous 
Section and is shown as dashed-dotted lines in Figure [B] 
The two formation models for MHs (ES and RE) represent 
two extreme effects of the IGM heating on the 21 cm signal 
from MHs. In each model, the signal of MHs does not vary 
as a function of IGM heating, so these lines are the same 
as the solid lines in Figure |4] STt from MHs is roughly pro- 
portional to the total mass of MHs, while it is not sensitive 
to the shape of their mass function. On the other hand, the 
IGM contribution strongly depends on its thermal state, and 
we present results for three different cases. 

First, we discuss the case in which only heating from 
ionizing radiation is present. Over the redshift range cov- 
ered by the simulations, the coUisional coupling in the mean- 
density IGM considered here is very weak. Thus, in the ab- 
sence of other decoupling mechanisms, the spin temperature 
of neutral hydrogen is very close to Tcmb and the 21 cm 
line is only slightly visible in absorption (dashed lines in 
Fig. [6|. However, as the flrst sources of radiation turn on 
and reionization begins, the energy of the ionizing photons 
heats the IGM in the vicinity of the sources to a significant 
high temperature (in this work we follow C06 and use the 
reference value lO^K, which is a reasonable a pproximation 
in regions ionized by stellar type s ources; see ICiardi et al.l 
l200ll : iPritchard fc FurlanettoH2007l ) . Thus the IGM in such 
regions may emit 21 cm radiation until it has been fully ion- 
ized, so that the mean STt changes from a negative to a pos- 
itive value. However, as reionization proceeds, the amount 
of neutral hydrogen decreases and the mean STt drops to 
zero again. 

In the extreme suppression case, we find that even in 
the presence of only coUisional decoupling and UV heating, 
the 21 cm signal from minihalos dominates the total signal 
only in the very early stages of reionization. This is because 
reionization suppresses the formation of new minihalos and 
photoevaporates existing ones, reducing their contribution 
with time. On the other hand, reionization increases the 
temperature of the partially-ionized IGM and the density of 
free electrons. In high temperature, partially ionized regions, 
H-e~ and H-p collisions are more effective than H-H colli- 
sions in coupling Ts to Tfc. As a consequence, the fraction 
of flux from partially ionized gas is higher than its volume 
filling factor, and the signal from the IGM increases with 
time until a volume averaged ionized fraction of xmi ~ 0.5 
has been reached, a.t z ~ 15, which corresponds to a peak of 
the emission. 

On the other hand, in the reformation model, the fiux 
from the IGM is reduced by the significant value of /co11,mh 
and pushed towards lower redshift by the delay in the reion- 
ization process. In this case, the period during which the fiux 
from minihalos is larger than that from the IGM is extended 
to Az ~ 4, and reaches a peak at z ~ 12. 

It should be stressed that the signal from the IGM dom- 
inates over that from the MHs because of the collisions be- 
tween atoms and free electrons and protons (we discuss the 
uncertainties of this in Sec. I3.3|l . In fact, if we ignore this 
contribution, the peak of STt from the IGM is only « 2 mK 
and the signal from MHs dominates down to z ~ 15.5 in 



The contribution of the IGM and minihalos to the 21 cm signal of reionization 7 




Figure 5. A map of 5Tt, in units of mK in the pure UV heating case and reformation case (see text for details). From top to bottom, 
the three rows correspond to the IGM, MHs, and the total. In each row, from left to right, the panels show slices at redshift 18.5, 14.6, 
12.8, 10.6 respectively. The comoving thickness of the slice corresponds to the dimension of one cell in the simulation, i.e. 223 kpc. 



the extreme suppression model, and at almost all redshifts 
in the reformation model. 

Next, we consider a case that includes the effects from 
a Lya background, leading to the mean STt evolution shown 
in Figure [7] Here we see that the Lya background produced 
by the sources in the simulation is strong enough to couple 
Ts to Tfc and to heat the IGM above the adiabatic temper- 
ature, but the mean temperature is still below Tcmb. This 
induces a strong mean signal in absorption, and a minimal 
emission is visible only at a late evolutionary stage in the 
reformation model. This is because reionization proceeds to 
lower redshift, when more sources form and produce a larger 
value of the Lya background. Yet even in this case the bulk 
of the emission is due to the partially-ionized regions. 

The reason for the very different behavior compared to 
the no Lya background case is the following. Without Lya 
photons the contribution to STt consists of a small negative 
signal from the neutral, cold IGM (because Ts « Tcmb), and 
a positive signal from the the partially ionized, hot IGM. 
In the presence of Lya photons instead, the dominant con- 
tribution comes from a strong negative signal from the the 
neutral, warm IGM (because Ts <C Tcmb) which exceeds the 
smaller positive signal from the partially ionized, hot IGM. 
In this case we find that the I GM signal always domin ates 
the total signal consistent with iFurlanetto fc OhI l|2006l ). 



Finally, we consider a situation in which Ts 2> Tcmb, 
resulting in the evolution of the mean STt shown in Fig- 
ure |8l As expected, in this case STb is always positive and 
it decreases as reionization proceeds and the available neu- 
tral hydrogen is depleted. Unlike the cases analyzed before, 
the signal corresponding to the extreme suppression model 
is always lower than that of the reformation model. This 
is due to the fact that, in this case, STi, is independent 
of Ts and depends only on the amount of neutral hydro- 
gen present at each redshift, which is higher in the refor- 
mation model. Moreover, the signal from the IGM always 
dominates over that from the MHs. This is because in both 
cases the condition Ts S> Tcmb is satisfied, meaning that 
(5Ti,,iGM oc /iGM.Hi and 5Tt,,MH oc /mh,hi, with a similar 
constant of proportionality. Thus /igm,hi 3> /mh,hi implies 

STi, JGM ^ STfjMH- 



3.2 rms fluctuations and pixel distribution 

We now calculate the rms fiuctuations in the brightness tem- 
perature distribution, {{STt — {STi,))^)^^^, as a function of the 
beamsize, AO, and redshift. We choose a slice corresponding 
to a bandwidth Av — 0.1 MHz, located at the center of the 
simulation box. A larger bandwidth (e.g. 1 MHz) would pro- 
duce lower fiuctuations because in that case inhomogeneities 
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Figure 6. Mean STi, evolution in the absence of heating sources, 
with the exception of the UV ionizing radiation. The solid, 
dashed, dashed-dotted and dotted lines refer respectively to the 
contribution of IGM+MHs, only IGM, only MHs and the volume 
averaged ionization fraction. The top panel refers to the extreme 
suppression model, while the bottom panel refers to reformation 
model. To guide the eye, we also plot a horizontal line correspond- 
ing to STi = 0. 



in the gas density and in the hydrogen ionized fractions 
would be more poorly resolved (see e.g. ICiardi fc MadaJ 

The results are shown in Figure l9l and 1101 for the three 
IGM thermal states and two MH formation models de- 
scribed in the previous Sections. Figure E] gives the rms as 
a function of beamsize at several different redshifts for both 
the reformation (upper panels) and the extreme suppression 
(bottom panels) model. Figure [TOl instead, gives the evolu- 
tion of rms for different beamsize in the same two models. In 
all cases, the signal increases as the angular scale decreases 
because the mass variance is larger on small scales. Another 
common feature is that, although the continuous formation 



Figure 7. As Fig. [6] but including heating by Lya photons. 

of minihalos increases the signal compared to the extreme 
suppression case, it does not increase the inhomogeneity of 
the gas, so the rms in the reformation model is similar to 
the rms of the extreme suppression model until the latter 
drops because of reionization. 

In the pure UV heating case, the signal from the IGM 
reaches a maximum at z ~ 14 (extreme suppression model) 
and at 2; 12 (reformation model), corresponding to the 
epoch when several neutral regions are still present, but the 
volume is roughly half ionized. Then it starts decreasing 
and drops dramatically when the IGM gets close to com- 
plete reionization. The signal from the MHs shows a simi- 
lar behavior, but in the reformation model the drop asso- 
ciated with photoevaporation is not evident because they 
are allowed to continuously reform. The general trend in 
the relative importance of the IGM and MHs contribution 
follows what already seen in the mean brightness temper- 
ature, i.e. the MHs contribution dominates in both models 
in the earliest stages of the evolution and when reionization 
is almost complete in the reformation model. Otherwise the 
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Figure 9. Expected rms brightness temperature fluctuations as a function of beamsize at different redshift. Upper and lower panels refer 
to the reformation and extreme suppression models, respectively. The three thermal states of the IGM are labeled in each panel. Thick 
(thin) lines refer to the IGM (MHs) contribution, at the redshifts 18.6, 14.6, 12.8, and 10.6. 




Figure 10. Expected rms brightness temperature fluctuations as a function of redshift with different beamsize. The meaning of each 
panel is the same as in Fig. [9] 



IGM contribution is dominant. In addition, the smaller the 
beamsize, the earlier the fluctuations that from IGM exceed 
those from the MHfl Thus, increasing the resolution of the 
beam might not help in detecting the signal from MHs (be- 
cause observations at higher redshift and lower frequencies 
are more difBcult), unless one chooses a region particularly 
rich in MHs. 



^ Here we only consider A0 = 0.1 because this corresponds to 
the size of one cell in the simulation. 



Also in the Ts ^ Tcmb case and in the Lya case, the 
behavior of the rms is similar to that of the mean brightness 
temperature, with the MH signal dominating only in the 
reformation model during the latest stages of reionization. 

We then calculate the pixel distribution of the bright- 
ness temperature, 5Ti,df /dSTt, as a function of STt for the 
IGM in the three cases and minihalos. Here df = Nsr^/Ntot, 
where Ntot is the total number of cells in the simulation box 
and NsT^ is the number of cells that lie in the bin dST},. In 
Figure [TT1 and [121 the panels correspond, from left to right 
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Figure 8. As Fig.|6]but in the Ts > Tcmb case. 



and from top to bottom, to pure UV heating case, Lya heat- 
ing case, Ts 3> Tcmb case and minihalos, respectively. 

The amphtude of the pixel distribution indicates the rel- 
ative contributions from different STi, bins to the total signal. 
The double peaks associated with the two-phase structure of 
the IGM (totally neutra l and partially ionized regions) are 
very prominent (see also ICiardi fc Madaxill2003l l . For exam- 
ple, in the Lya heating case, at an early stage (dashed lines), 
when most of the signal comes from completely neutral re- 
gions, only one peak is visible in absorption, corresponding 
to 5r6,peafc ~ —150 mK. However, as ionization proceeds, the 
amplitude of the peak decreases and moves toward a positive 
value; meanwhile, a new, smaller peak appears with positive 
value at STi, ~ several mK, which comes from the partially 
ionized regions (see the dashed and dashed-dotted lines in 
Fig, inl and Fig. I12p. The three redshifts chosen show three 
typical situations: when the totally neutral regions domi- 
nate, when the partially ionized regions begin to produce a 
significant signal and when the contribution from partially 
ionized regions is more evident, respectively. The behavior of 
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Figure 11. Pixel distribution for both the IGM and minihalos 
for the extreme suppression model. Different lines correspond to 
different redshifts. In each IGM panel, solid, dashed and dashed- 
dotted lines mean redshift 18.1, 15.5, 14.8 respectively. In the 
minihalos panel instead, pixel distribution is given at redshift 
16.6 (solid) and 13.7 (dashed). 



the pixel distribution is very similar for the extreme suppres- 
sion and the reformation model, the only difference being a 
delay in the reformation model. In the S> Tcmb case, 
the trend is the same as in the Lya heating case, i.e. as 
ionization proceeds, the peak corresponding to the totally 
neutral regions decreases and moves toward zero, while the 
peak originating from the partially ionized regions increases. 
However, the situation is slightly different in the pure UV 
heating case, in which the contributions from the partially 
ionized regions become significant at an earlier stage. 

Finally, the dominant minihalo contribution at all red- 
shifts comes from cells with STb « 1-2 mK, which is a typical 
value of beam-averaged differential brightness temperature 
in our simulation (this value depends on the c osmological 
param eters and the matter power spectrum, see llliev et al.l 
l2002h . Unlike the IGM signal, here there is always only one 
peak; however, since we can measure only a total flux, this 
feature is not useful to distinguish the signal of minihalos 
from that of the IGM. The sawtooth features come from the 
discreteness effect introduced when calculating /mh, which 
shows also in the spikes of the IGM signal. 

3.3 Model uncertainties 

In this Section we discuss the main approximations and un- 
certainties of our model and their impact on the results, in 
addition to those already discussed in C06. 

• In the pure UV heating case, partially ionized regions 
contribute significantly to the emission signal. However, we 
note that the resolution of our simulation is only « 223 kpc, 
which is not enough to resolve the I- front. This leads to the 
possibility of overestimating the signal from partially ion- 
ized regions by mixing together totally ionized and totally 
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Figure 12. Same as Fig. 1111 but for the reformation model. In 
each IGM panel, solid, dashed and dashed-dotted lines mean red- 
shift 18.1, 14.2 and 12.2 respectively, while in the minihalos panel, 
solid line refers to redshift 16.6 and dashed line refers to redshift 
11.1. 



neutral regions. In fact there are many uncertainties related 
to the sub-grid neutral filling factor, including the inhomo- 
geneity of the densi ty field on small scales and the spectrum 
of the sources (e.g. Illiev et al.|[2006l ). 

In our work we refer to "partially ionized" cells as regions 
heated by UV photons that still contain enough neutral hy- 
drogen to emit a 21 cm signal. Note that this somewhat 
arbitrary definition is used only to estimate the relative con- 
tribution of gas in different physical conditions, while all the 
cells in the simulation have been used to calculate the total 
21 cm signal. In all our previous calculations, we identi- 
fied partially ionized cells as those with ionization fraction 
X above a threshold value of 0.1, but a value of 0.01 leads 
to similar results. If x is very close to 1.0, the signal from a 
cell is negligible and the cell is referred to as totally ionized. 

If we assume that in a partially ionized cell only a fraction 
fpi of the gas is partially ionized with an ionization fraction 
of Xpi, while the rest is either neutral and cold, or ionized 
and hot and thus irrelevant for the contribution from the 
partially ionized regions. Furthermore, we should consider 
that only a fraction of gas fpi x (1 — Xpi) is contributing 
to the signal. Both fpi and unknown, and vary in a 

large range cell by cell. As a reference, if we assume Xpi — x 
we find that at the redshift when half of the gas has been 
ionized, the mean STt of the total simulation volume for 
fpi equal to (0.0, 0.001, 0.01, 0.05, 1.0) are (-0.054, -0.051, 
-0.012, 0.21, 5.63) mK respectively for the extreme suppres- 
sion case. Note that {fpi, STt) = (1.0, 5.63 mK) corresponds 
to the peak of the dashed line in upper panel of Fig. |6] 

So if the sub-grid neutral filling factor is much lower than 
1.0, we expect no emission mean signal or non-dominant 
partially ionized IGM signal in the pure UV heating case. 
However, the pure UV heating model ignores the pre-heating 
by X-rays and Lya photons, thus we expect it to be appro- 



priate to describe only the early stages of reionization. Also 
the small emission in the Lya coupling case suffers from this 
uncertainty, although our main trends and conclusions will 
not change. Finally, in the strong heating case, the infiuence 
of resolution should be negligible. This is natural since in 
those two models, the dominant signal is from the neutral 
gas. 

• To derive the mass evolution of MHs undergoing photo- 
evaporation we use the fitting formula bv Illiev et al.l l|2005h . 
The formula reproduces the results from their numerical 
simulation, but it is approximate. Moreover, it is accurate 
only for a constant ionization fiux, while in the simulations 
by C06, the fiux assumes different values at different times 
and locations. We tried two different methods to derive 
eq. (Ill|l in the presence of a time dependent flux, correspond- 
ing to eq. (|12|l and eq. (|14|) . which consistently predict the 
5Tb within few percent. As in C06, for most halos, the evap- 
oration time is comparable or smaller than the time between 
two outputs (for example, in the extreme suppression model 
~ 66% of the cells lose 80% of their MH mass between two 
outputs). We also do not expect the use of those formulae to 
introduce uncertainties in addition to those already present 
in the C06 implementation. 

• To calculate the 21 cm flux from an ind ividual mini- 
halo we use the formula from Illiev eTall (|2002l ') and then we 
integrate over the mass function in eq. (I10|l to obtain the 
beam-averaged effective differential brightness temperature 
of eq. (jSj. This formula is only suitable for halos that are 
in post-collapse equilibrium with a TIS profile. Yet, when 
a halo is being photoevaporated, both its density and tem- 
perature profile may be substantially differenllj. The new 
profile should depend on the local ionization flux. If we de- 
fine the dynamical timescale as t^yn ~ Trrris/vc, where Vc 
is the halo circular velocity, for a halo with mass 10^ Mo 
at redshift 15, tdyn ~ 150 Myr while t^v ~ 230 Myr ii 
Fo = 1. So the timescales are comparable and the halo pro- 
file could be chan ged during photoevaporation. In fact, also 
Illiev et all (|2005l ) notice that the typical evaporation time 
is long enough to see a gas-dynamic back-reaction. Since 
the effects of this change on the 21 cm signal of MHs is 
still not clear, we ignore the effect in this work. In addition, 
heating of the IGM can also change the profile of a halo, 
mak e it less concentrate d and thus increase the 21 cm sig- 
nal (|Oh fc Haimanll2003l ). On the other hand, the equation 
Ts 2> TcMB will always be satisfied, so the impact of the 
temperature change on the 21 cm signal should be much 
less than that of the change on the density profile. 

4 CONCLUSIONS 

We have studied the 21 cm signal from both the IGM 
and minihalos base d on the reionization simulations run by 
ICiardi et al.l (|2006l ). These include a "sub-grid" prescription 
to take into account the photoevaporation of MHs in two 
different models for their formation (extreme suppression 
and reformation). Our estimates follow the changes in the 
mass function of MHs and consider three heating case of the 
IGM. We find that: 

And in general, also if undisturbed, the halo c ould follow a 
different shape (e.g. iNavarro. Frenk fc Whitelll997t) . 
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• the 21 cm signal tightly relates to the thermal history of 
the IGM. In all the three configurations (heating purely by 
the UV photons, Lya and UV heating, and strong heating), 
the IGM dominates the 21 cm signal at almost all redshifts, 
except at a very early stage when the gas is still cold and a 
Lya background has not been constructed or at the end of 
reionization in the reformation model, when almost all the 
remaining neutral gas is in MHs. The small contribution of 
minihalos to the total 21 cm flux is partially due to the fact 
that the collapsed fraction of MHs is relatively small at high 
redshift. Minihalos are signiflcant emission sources in the 
absence of IGM heating, however, in partially ionized regions 
which are heated by UV photons, where there is enough 
residual neutral gas and the temperature is much larger than 
TcMB , the IGM also provides a contribution to the emission 
signal. In the presence of a Lya background, neutral gas 
will be visible in strong absorption until a later stage when 
neutral regions are rare. In such case, totally neutral regions 
contribute more than partially ionized regions and MHs to 
the total 21 cm signal. 

• Although MHs are directly visible in 21 cm signal only 
under particular circumstances, they nevertheless play a sig- 
nificant role in affecting the ionization state of the IGM 
and the corresponding 21 cm flux. Their presence can delay 
reionization by as much as Az « 4, depending on the effect 
of feedback on their formation (COG). More speciflcally, in 
the pure UV heating case, the 21 cm signal from the IGM 
in the reformation model is weaker than that in extreme 
suppression model at the same redshift, since MHs delay 
the photoheating associated with reionization. In the strong 
heating case on the other hand, the 21 cm emission is higher 
in the reformation model than in the extreme suppression 
model at the same redshift, because absorption from MHs 
leaves more neutral gas in the IGM. This difference shows 
that in the pure UV heating case, the temperature is the 
most important factor which can impact the 5Ti, of the IGM, 
while in the strong heating case the most important factor 
is uhi- 

• While the 21 cm signal from MHs exhibits differences 
from the IGM signal, including the pixel distribution, de- 
tecting such features will not be within the capabilities of 
the planned generation of radio telescopes (e.g. LOFAE0, 
MWA0, 21cmi^. In principle, MHs are more abundant in 
rich cluster regions, which implies a higher 21 cm signal. In 
practice, in these regions also contain more ionizing sources 
that suppress MHs formation and evolution, reducing their 
signal. 

Conclusions similar to those derived from our simu- 
latio ns have been reached by previous analyt ic estimates 
(e.g. 'Oh & Mack''2003'; 'Fur lanetto fc Ohl l2006h. For exam- 
ple, iPuHajiMto & Oh (2006) calculated the signal from both 
the IGM and the minihalos, using a semi-analytic approach 
that includes the X-ray heating and feedback effects on halo 
formation. They found that in the presence of Wouthuysen- 
Field coupling, MHs are always buried by the IGM signal, 
even if the Lya bac kground is small. On the other hand, 
IShapiro et al] l|2006l) concluded that before reionization but 

* http://www.lofar.org/ 

^ http:/ /www. haystack. mit.edu/mwa/index. html 
^ http://21cma.bao.ac.cn/ 



after z ~ 20, ignoring the radiative processes but taking 
into account shock heating, minihalos could dominate the 
total signal. This conclusion is consistent with our pure UV 
heating case at the early stage. 
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